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We study heat and charge transport through a normal diffusive wire coupled with a supercon- 
ducting wire over the region smaller than the coherence length. Due to partial Andreev reflection of 
quasiparticles from the interface, the subgap thermal flow is essentially suppressed and approaches 
zero along with energy, which is specific for diffusive structures. Whereas the electric conductance 
shows conventional reentrance effect, the thermal conductance k rapidly decreases with tempera- 
ture which qualitatively explains the results of recent experiments. In the Andreev interferometer 
geometry, k experiences full-scale oscillations with the order parameter phase difference. 
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Manifestations of the proximity effect in the electron 
transport of the hybrid normal metal-superconductor 
(NS) structures are in the focus of current extensive re- 
search. Up to now the electric conductance in NS-hybrids 
has been receiving much more attention than thermal 
conductance; the imbalance can be ascribed to difficul- 
ties one encounters in carrying out thermal transport ex- 
periments in mesoscopic samples. A remarkable break- 
through in recent measurements 0, 0] of both the ther- 
mal conductance and thermopower in an Au (N) diffusive 
wire of the micron length and submicron cross-size and 
thickness, where an Al (S) needle-like sample of similar 
parameters is deposited across it, as shown in Fig. ^a, 
calls for and motivates a detailed theoretical investigation 
of heat and charge transport in mesoscopic NS-hybrids. 
In our Letter we develop a theory of thermal and elec- 
tric conductances in mesoscopic proximity structures in 
the diffusive limit, which, as we show below, describes 
the experimental situation. We find that while the elec- 
tric conductance exhibits the conventional reentrance be- 
havior, the thermal conductance rapidly decreases with 
temperature, as k ~ T 4 , in a qualitative agreement with 
experimental findings. 

Experiments 0, 0| showed a tiny (within a few per- 
cents) change in the electric conductance of Au due to 
the proximity effect in accordance with the past studies. 
The thermal conductance, in a contrast, dropped with 
temperature decreasing, by the order of magnitude as 
compared to its value kn(T) in the normal state. This 
result seems to excellently follow the original idea by An- 
dreev Q first applied to the thermal conductance of the 
intermediate state of a superconductor Q: If the inter- 
face resistance is negligibly small, and the wire thick- 
ness is smaller than the coherence length then the 
sandwich-like NS contact in Fig. ^ a can be modelled b 
an inset of the superconductor into the N-wire, Fig.Q]b 
Then the superconductor lead in Fig.^b plays the role of 
a quantum barrier for normal quasiparticles with the en- 
ergies E smaller than the superconducting order param- 
eter A. Low-energy electrons hit the barrier and convert 
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FIG. 1: Two models of the NS interface: sandwich geometry 
(a), and inset geometry (b). 
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FIG. 2: (a) Upper view of the model of the NS interface in 
Fig. with the ballistic trajectories of the electrons (solid 
lines) and holes (dashed lines) ; (b) the geometry of house-like 
Andreev interferometer with the magnetic flux 

into the retro-reflected holes (Andreev reflection), carry- 
ing the same energy back in the opposite direction. Thus, 
the subgap thermal flow through the superconductor is 
blocked, and only quasiparticles with energies E > A 
participate in the heat transport. At low temperatures, 
T -C A, the thermal conductance follows nearly expo- 
nential temperature dependence, k ~ Texp(— A/T) 
used in Ref. for data fitting. 

Given a very satisfactory fit, the above line of rea- 
soning would have been conclusive, unless the disturbing 
observation that in the experiments of 2] , the cross size 
w of the N- and S-leads was considerably smaller than 
£oi making the Andreev barrier semi-transparent for the 
subgap particles. This raises serious doubts about the 
applicability of the above 'macroscopic' consideration, re- 
quiring w ^> £oi t° the quantum problem, w < £0, where 
the subgap quasiparticles may well dominate the thermal 
conductance at low temperatures. This is the case in bal- 
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listic structures, with the elastic scattering length I much 
larger than £o and w , where the heat flow is carried by 
quasiparticles traversing the proximity region and over- 
coming the weak Andreev barrier (see Fig.|2a). Such ef- 
fect was considered within a model of a small S-inclusion 
in the ballistic N-wire Q. In the geometry of Fig. [21a, the 
quasiparticles which propagate into the S-lead over the 
distance £o, undergo complete Andreev reflection, and 
sink out from the heat transport 0. In the limit of 
almost transparent Andreev barrier, w <C £o> the ratio 
k/kn in the ballistic structure is basically determined 
by the interface geometry and therefore should be tem- 
perature independent, except in a vicinity of T c , where, 
as £o becomes larger than the length of the S-lead, the 
Andreev reflection is suppressed, and k returns to kjv- 
In this Letter, we focus on the diffusive limit, £ <C 
that describes the experimental situation in Refs. 
and reveals a wealth of interesting physics. In this 
case, due to multiple coherent backscattering of elec- 
trons by the impurities within the proximity region Q, 
the number of Andreev reflected quasiparticles effectively 
enhances. This enhancement becomes especially pro- 
nounced as the energy approaches zero, and, correspond- 
ingly, the size £n (E) = yjKD/2E of the proximity region 
in a normal metal with the diffusion coefficient D in- 
finitely increases. This effect is known as the reason for 
the zero-bias conductance peak in NS structures with an 
opaque interface ||. In our context, enhanced Andreev 
reflection means the suppression of the heat flow at small 
energies even for the weak Andreev barrier, w <C £o> and, 
as we show below, results in a rapid decrease of k/kn 
with temperature as observed in the experiment. Fur- 
thermore, the temperature dependence of k appears to 
be power-like rather than the exponential, in contrast to 
the model of the "impenetrable" Andreev barrier 0, Q 
which blocks the heat transport within the entire sub- 
gap region. Due to high sensitivity of the heat flow to 
the proximity effect, one can also expect essential depen- 
dence of the thermal conductance on the order parame- 
ter phase difference <p applied to the proximity structure 
(in the ballistic case, such dependence was predicted in 
Ref. 0). A hallmark of the phase-coherent heat transport 
was found in thermopower measurements 0, E| . 

We consider both electric and thermal conductances 
of the proximity structure in Fig. [21a in a diffusive limit, 
assuming the N-lead of length Ln {—Ln/2 < x < Ljv/2) 
to be connected to normal reservoirs, and the S-lead of 
length Ls (0 < y < L$) attached to the middle of the 
N-lead. In such geometry, the quasiparticle flow along 
the S-lead is blocked, and therefore the quasiparticle dis- 
tribution in this lead is spatially uniform. For this rea- 
son, it is enough to solve one-dimensional diffusive kinetic 
equations [3, LLiJ for the distribution functions f± in the 
N-lead, d x (D±d x f±) — 0, neglecting spatial variations 
across the leads at w <C £o- The diffusion coefficients D± 
are defined through the retarded branch of the spectral 



angle 9 as D + — cos 2 Im(9 and Z?_ = cosh 2 Re6>. The 
spectral angles 9m and 9s in the N- and S-lead, respec- 
tively, are to be determined from the Usadel equation 



2{E sinh 9 - A cosh 9) = ihVd 2 9, 



(1) 



with the boundary conditions 9n — at the normal reser- 
voirs, and d y 9s = at the edge of the S-lead. At the 
transparent NS interface, the spectral angle is continu- 
ous, 9n(0) = 9s{0) = #o, and obeys the current conserva- 
tion law [13, gsd y 9 s (0) = ±2g N d x 9 N (T0), where g N . s 
are the conductances of the leads per unit length. 
The kinetic equations have the first integral, 



D±d x f ± = I ± (E), 



(2) 



where the spatial constants I±, which have the mean- 
ing of the spectral densities of the probability and elec- 
tric currents, respectively, are to be determined from the 
boundary conditions for f±. The electric current / and 
the thermal flow Q in the N-lead are related to I± as 

/>oo />oo 

I=2K I_{E)dE, Q = ^ EI+(E)dE. (3) 
e Jo z Jo 

The thermal flow arises at different temperatures T\ 
and Ti of the N-reservoirs, which imposes the bound- 
ary conditions for the quasiparticle density function 
/ + (=FLjv/2) = tanh(_E/2T 12 ), where a small thermo- 
power is neglected. Assuming the temperature difference 
to be small, T x - T 2 < T = (T x + T 2 )/2, Eqs. © and © 
give the expression for the thermal conductance, 
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E 2 Rq(E) dE 
Ln{T\ - T 2 ) ~ 2tt 2 T 3 J q cosh 2 (^/2T) ' 



(4) 



where Rq(E) = (D+ l {E,x)) is the dimensionless spec- 
tral thermal resistance, and the angle brackets denote 
average over the length of the N-lead. 

In the problem of charge transport, we assume the left 
and right reservoirs to be voltage biased at =F^/2 and 
maintained at equal temperature T, which results in the 
boundary condition for the charge imbalance function /_ , 
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(5) 



Using Eqs. (J2J), © and ©, we obtain the zero-bias con- 
ductance of the N-wire, 
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Rj 1 (E)dE 
cosh 2 (£/2T)' 



(6) 



where Gn — 9n/Ln is the normal conductance, and 
Ri{E) = (DZ 1 (E,x)) is the spectral electric resistance. 

The calculation of the spectral angle 9n(E,x), which 
enters the coefficients D± , is to be performed numerically. 
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However, neglecting spatial variations in A within the 
proximity region (that may result in some non-crucial 
numerical factors at most) and assuming the lengths of 
the leads to be much larger than the characteristic scales 
£/v S of spatial variation of the spectral angle, 
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(7) 



one can apply the following analytical solutions of Eq. 
for the structure with semi-infinite leads, 



tanhM^ =taim ^) e -H/^ 

A J 



9 = \n 




e -y/is 
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(8) 
(9) 



L>N \9S J L>S 

Here 9 b = Arctanh ( A/E) is the spectral angle in a 
bulk superconductor, and the parameter r determines 
the strength of the proximity effect (see comments to 
Fig. |2J. Similar solution for the normal part of a long 
SNS structure was used in Ref. 

The semi-infinite-lead approximation enables us to es- 
timate the magnitude of k in the most interesting case 
of small temperatures, T <C A. Although the coeffi- 
cient D + is suppressed within the whole proximity region, 
M ^ Cat ~ £o\/A/T, the contribution 5Rq of this region 
to the thermal resistance Rq comes from the narrow (of 
the order of £o = ^KDm/2A) vicinity of the crossing 
point, in a contrast to the proximity correction to the 
zero-bias electric conductance which is formed over the 
much larger scale £/v. Within this region, the coefficient 
D + is anomalously small and turns to zero at E — * 0, 



(E/A)(\x\/Co+VP)'' 



(11) 



which results in the following estimate of 5Rq at the 
characteristic energies E ~ T, 



6R Q (T)^T /T, T = A£o/L 



JV\ 



(12) 



Thus, at T < To, the proximity region dominates the net 
thermal resistance giving rise to the power-like decrease 
of k ~ k n T/T ~ T 4 with T. At very low temperatures, 
smaller than the Thouless energy E^h = KDn/L 2 n , the 
N-reservoirs begin to affect the quasiparticle spectrum at 
the contact area, and the approximation of semi-infinite 
N-lead fails. In this case, -&rh provides the cutoff in the 
decrease of the coefficient D + at small E, which therefore 
saturates at D + ~ Et^/A. As the result, the thermal 
conductance at T < i?Th starts to decrease slower, as 
kn ~ T 3 , with a small prefactor Et^/Tq ~ ^o^/t/Ln- At 
the same time, the electric conductance shows behavior 
typical for the proximity structures: as T decreases, G(T) 
approaches maximum and then returns to Gn- 
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FIG. 3: Temperature dependencies of k/km (a) and G/Gn 
(b) obtained by numerical solution of Eq. Q, in comparison 
with the Andreev model (dashed curves), at Lm — 10£o- 

Shown in Fig. |3| are the results of calculation of ther- 
mal and electric conductances within the whole relevant 
temperature region < T < T c for different magnitudes 
of the parameter r. If the conductivity of the N-lead 
is small, r C 1, the spectrum within the S-lead is al- 
most unperturbed [9q « 9b, see Eq. IjlUI) ]. and the cross- 
ing area plays the role of a superconducting reservoir. 
In terms of quasiparticle motion, quasiparticles 'would 
prefer' to diffuse into the high-conductive S-lead, where 
most of them undergo Andreev reflection thus addition- 
ally reducing k and increasing G. In the opposite limit, 
r» 1, the quasiparticles avoid penetration into the low- 
conductive S-lead, which suppresses the proximity effect 
(#o ~~ * 0), and, correspondingly, enhances k and reduces 
the peak in G(T). Note that for all reasonable values of 
r, our quantum model gives the magnitude of k consid- 
erably larger than that due to the classic model 0, HI • 

The high sensitivity of thermal conductance to the de- 
tails of the quasiparticle spectrum can be used for a study 
of the effects of phase coherence in the mesoscopic devices 
known as Andreev interferometers. Usually, the interfer- 
ometer circuits contain a superconducting loop with the 
order parameter phase difference <f> controlled, e.g., by 
the magnetic flux, and connected to an SNS junction. 
The electric conductance of the junction (or the N-wire 
attached to the junction) depends periodically on cf>. 

The basic features of the phase-coherent heat trans- 
port in the proximity structures can be demonstrated 
within the model of the "house interferometer" 0, Q 
shown in Fig. [21b. Assuming the length of the junction 
arms to be smaller than £o, and the normal conductivity 
of the S-loop to be much larger than the conductivity of 
the N-wire, one can neglect the change in the junction 
quasiparticle spectrum due to proximity to the N-wire. 
In this limit, the spectrum at the contact area is simi- 
lar to the BCS-like spectrum of a separate short diffu- 
sive SNS junction, with the phase-dependent energy gap 
E g (<j>) = A|cos(0/2)| 0. This implies that the con- 
tact area behaves as an impenetrable Andreev barrier 
with the phase-dependent height E g (</>), and therefore 
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FIG. 4: Phase dependencies of the thermal (a) and electric 
(b) conductances in the geometry of house interferometer in 
FigEDb, for different temperatures (Ljv = 10£o)- 

the heat transfer is governed by the quasiparticles with 
the energies E > E g {<j)). Within this region, we apply 
Eq. J3J, thus taking into account spatial variations in 
the diffusion coefficient D + within the proximity region. 
These variations are a diffusive analog |LJ] of the bal- 
listic over-the-barrier Andreev reflection from a rapidly 
varying order parameter potential. 

The phase dependence of the thermal conductance is 
presented in Fig.^Ja. Due to Andreev reflection from the 
contact area, k is generally suppressed, according to the 
mechanism discussed above. However, as the phase dif- 
ference approaches odd multiples of tt, the energy gap E g 
closes, and all quasiparticles freely diffuse through the N- 
wire, which restores the normal value of n. Thus, at low 
temperatures, k(4>) exhibits strong oscillations with the 
phase </>, with sharp peaks at 0mod27r = tt. The elec- 
tric conductance shows oscillations with sharp dips at the 
same points, but with much lower amplitude (Fig. 0Jb). 

In conclusion, we have developed a theory of the heat 
transport through the normal part of diffusive NS na- 
nostructures. The Andreev reflection from the NS in- 
terface blocks the quasiparticle probability current, and 
therefore the heat flow through the normal wire, cou- 
pled with the superconductor over a small area, is es- 
sentially suppressed. Our approach takes into account 
self-consistently the influence of both the normal metal 
and superconductor on the quasiparticle spectrum within 
the proximity region, extending thus the ballistic concept 
of partial Andreev reflection from the superconductor of 
finite width to diffusive proximity systems. The impor- 
tant feature of the diffusive structures is the essential en- 
hancement of the Andreev reflection probability at low 
energies due to multiple returns of coherent quasiparti- 
cles to the contact area. For this reason, the proximity 
region dominates the net thermal resistance of a long 
normal wire with the length ijy 3> Co at l° w tempera- 
tures, even if the size of the contact area is much smaller 
than the coherence length £q- This results in a power-law 
decrease in the thermal conductance with the tempera- 
ture, k *~ T 4 , transforming into a cubic law at the tem- 



peratures smaller than the Thouless energy. The effect 
becomes more pronounced when the carriers mobility in 
the superconductor is higher than in the normal metal. 
If the wire is attached to the SNS junction with the order 
parameter phase difference <j> between the electrodes, the 
thermal conductance reveals full-scale oscillations with <fi 
and shows large peaks at odd multiples of tt. 

In the experiments Q , the results of measurements of 
k(T) were fitted by the Andreev formula |jg, which as- 
sumes complete blockade of the thermal flow within the 
entire subgap region, E < A, and therefore leads to the 
exponential temperature dependence of k. Although this 
fitting looks satisfactory, the physical background for the 
applicability of such a simple model remains unclear, be- 
cause of the small width w » 0.4£ of the wires and un- 
avoidable suppression of electron-hole correlations in the 
S-wire within the proximity region. This calls for fur- 
ther investigations, including independent measurements 
of k(T) in the normal state and the extension of theoret- 
ical calculations for rather complicated geometry of real 
proximity structure in Ref. 0- 
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